Study of Nondiagonal Parton Distribution Models 



Andreas Freund, Vadim Guzey 
Department of Physics, The Pennsylvania State University 
University Park, PA 16802, U.S.A. 



0^ 

Q> '. Abstract 

\ In this paper we examine predictions from different models of nondiagonal 

parton distributions. This will be achieved by examining whether certain 
predictions of relationships between diagonal and nondiagonal parton distri- 
J> butions also hold after having evolved the different distributions. 
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I. INTRODUCTION 

In recent years deeply virtual Compton scattering (DVCS), hard exclusive electropro- 
duction processes and the distributions which contain the nonperturbative information for 
those processes, namely the nondiagonal or nonforward parton distributions, together with 
their evolution have been a target of intense studies [p]-[21[ . 



Nondiagonal parton distributions, in particular, have attracted a lot of attention and 
numerical studies as to their behavior under their leading order (LO) evolution were con- 
ducted recently 0-0. Also recently, the next-to-leading order (NLO) generalized Efremov- 
Radyushkin-Brodsky-Lepage anomalous dimensions were computed employing conformal 
symmetry arguments [|I^ making NLO studies in the near future possible. The NLO kernels 
for the DGLAP region (xi > A) of the nondiagonal parton distributions can be reconstructed 



via the method described in Ref. 22 



In the above mentioned numerical studies, a particular ansatz for the initial distribution 
for nondiagonal evolution was made, namely that the nondiagonal and diagonal initial dis- 



tributions were taken to be equal in the normalization point. As was pointed out in ||23 
this ansatz cannot be justified except at very low normalization point Qo where the parton 
distributions are essentially flat. In the following, we will explore a different ansatz and 
investigate whether predictions based on this ansatz hold after evolution. It is crucial for 
the predictive power of theoretical calculations of exclusive, hard diffractive processes to 
have sensible nonperturbative models of nondiagonal parton distributions as long as they 
have not been extracted from experiment. 

The paper is structured the following way. In Sec. H we recap a few basic definitions 



about nondiagonal parton distributions and the involved kinematics, in Sec. Ill we will briefly 



explain the model under consideration as derived by Radyushkin ||23| and the predictions 
about the relationship between nondiagonal and diagonal parton distributions, in Sec. ^ we 
will give the results of our study, making some remarks about the behavior of nondiagonal 
parton distributions in NLO and, finally, in Sec. |V| we give concluding remarks. 



II. BASICS 

The kinematics as encountered in the appearance of nondiagonal parton distributions 
can be most easily described by the DVCS process 7*(g) + P{p) iW) + P{p')- The fact 
that the initial and final state protons have different momenta leads to the introduction of 
nondiagonal parton distributions. The nondiagonal character of these distributions can be 
seen in the definitions of nondiagonal parton distributions as Fourier transforms of matrix 
elements of bifocal, renormalized light-cone operators (see Ref. and below) , the light- 
cone operators being sandwiched between states of different momenta as compared to the 
diagonal case where the states have the same momentum. 

The important kinematical parameters are the following: is the virtuality of the 
probing photon, t = {p — p'Y is the momentum transfer to the final state proton, XBj is the 
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P(p) P(p') 

FIG. 1. The lowest order handbag contribution to DVCS with = —q^ and q'"^ = 0. 

usual Bjorken scaling variable, A = xi — X2 is the asymmetry parameter of the process and 
Xi and X2 are the longitudinal light-cone momentum fractions of the partons connecting the 
intermediate state to the hard scattering part (see Fig. 1). Note that in the case of DVCS 
XBj = A. 

For our study we need the definitions of nondiagonal parton distributions in terms of 
light-cone operators, their evolution equations and LO evolution kernels which we will review 
next. 

The definitions of nondiagonal quark and gluon distributions are given by 
dy- 
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with the following evolution equations for the singlet (S) and nonsignlet (NS) case 
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with the generahzed + regularization prescription necessary for the nondiagonal case which 
is described in detail in Ref. @,|^J^. 



III. RADYUSHKIN'S MODEL 



In the following we will briefly review the simplest model for an initial nondiagonal parton 
distribution as proposed by Radyushkin in Ref. |2^. The starting point are double distri- 



butions F{x,y) = F{x,y,t = 0) which behave like a distribution amplitude in the variable 
y and like a parton distribution in the variable x. They are expressed through multivariable 
Fourier transforms of matrix elements of bilocal, renormalized, light-cone operators [|2|,|23|. 
Symmetry and spectral conditions suggest the following ansatz for the double distribution 



such that with the normalization condition 

rl—x 

h{x) = / h{x,y)dy , (5) 
Jo 

one obtains the diagonal distribution /(x) 

f{x)= r^F{x,y)dy. (6) 

•J 

The simplest realization of the ansatz Eq. (^) is 
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F{x,y,Qo) = S{y-k{l-x))f{x) . (7) 

The 5-function in y prevents any spread of the double distribution in the y direction and k 
is a number smaller than 1 and is determined by the powerlaw behavior of parton densities 
at small x. 

In order to obtain the nondiagonal parton distribution in the DGLAP region {xi > A) 
from the double distribution, one has to integrate over y with an explicit A dependence 

A, Qo)i = I Fi{xi - A|/, y, Qo)dy . (8) 
Jo 

Using Eq. (|^) one finds the following initial distribution which will be our input for the 
evolution 

F.(.,.A.O„) = ^/,(^.0„) , (9) 



with i = q,q, g and we chose k = 0.5 which is in line with the arguments in ||23| that the crest 
of the double distribution is shifted towards the line y = {1 — x)/2. Note that the ansatz 
of the previous numerical studies |]^,|^ was that the diagonal and nondiagonal distributions 
were equal in the normalization point which corresponds to a double distribution similar to 
Eq. ® with 6{y - k{l - x)) 6{y). 

The above made ansatz leads to the following predictions (see Ref. for more details) 
for nondiagonal parton distributions: 

P(n)- ^^^^'^'^^ ^, (A/2)/,(A/2,g) 

""^^^ A/,(A,Q) - A/,(A,g) ^'^^ 

and 

A, Q) ^ {x, - A/2)f,{x - A/2, Q) . (11) 

The first equation is of importance in the case of DVCS since there one has xi = A = xbj- 
The second equation supposedly holds for any Q. If the above predictions hold true even 
after evolution it would give a very useful approximation of nondiagonal parton densities. 
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IV. RESULTS OF EVOLUTION 



Our input distributions for the diagonal parton densities f{x) in Eq. (0) will be CTEQ4M 
and CTEQ4LQ The reason why we chose these distributions is simple. Both CTEQ4M 
at a fairly high normalization point of Qo = l-Q GeV and CTEQ4LQ at a very low normal- 
ization point of Qq = 0.7 GeV are discernably different in the initial shape. This will help 
us discriminate fairly easily how robust to the initial shape the predictions of the previous 
section are. 

Our input distributions for the nondiagonal evolution are those satisfying relation (^. 
Since the nondiagonal evolution kernels are known explicitly to leading order (LO) only, we 
only use LO kernels for the diagonal evolution also, in order to have a consistent comparison 
between the nondiagonal and diagonal case. 

Let us first discuss how well Eq. (|Tl|) is satisfied. The results of our numerical study are 
given in Figs. (^,^,^,^). Although the input distribution for the nondiagonal evolution was 
chosen according to Eq. (|ll]) at Q = Qo, the evolution hardly violates this relationship at Q > 
Qq. One can see from Figs. that as the nondiagonal density evolves with Q the ratio of 
the nondiagonal to diagonal gluon parton density g{xi, A, Q) /{{xi — A/2)G'(xi — A/2, Q)) 
stays within 10% for all Xi and Q, for both CTEQ4M and CTEQ4LQ. Note that in these 
figures Ximin = I.OOOIA. As Xi increases the difference between the nondiagonal and diago- 
nal densities becomes small. This is a natural behaviour of the nondiagonal densities since 
for xi ^ A all the asymmetry related effects are unimportant. 

Next, we present the ratio of the nondiagonal to diagonal quark parton densities 
q{xi, A, Q) /{{xi — /S./2)Q{xi — A/2, Q)) in Figs. (^,||). The quark parton density is defined 

Q(x, Q) = u{x, Q) + Q) + s{x, Q) + u{x, Q) + d{x, Q) + Q) . (12) 

Here we observe a similar tendency as in the case of gluons - at xi close to A - the difference 
between the nondiagonal and diagonal densities is different from 1. However, in the case 
of quarks the ratio is significantly larger than 1. In fact, the nondiagonal quark is about 
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3 times larger than the diagonal one at xi ~ A (here ximin = 1.0001 A) [[ This result is 
not too surprising, in light of the findings of Ref . |^ , which showed a large deviation of the 
nondiagonal quark distribution from the diagonal one for a;i ~ A with a strong enhancement 
in the deviation for a low normalization point. As xi becomes significantly large than A the 
ratio quickly and smoothly approaches 1, as expected. 

To summarize this set of figures, we conclude that the prescription of Eq. (pUD , where 
one shifts the argument of the diagonal parton density by A/2, decreases the percentage 
deviation of the nondiagonal to diagonal parton density by approximately a factor of 3 — 4 
for gluons (compare to the relevant figures from Ref. [0]), giving a very good agreement for 
all X and Q. For quarks the approximation of Eq. ([TT|) is much worse as compared to the 
gluon case, however it becomes relatively good for xi > 2 to 3 A. 

Next we find that Eq. ([10|) is fulfilled with an accuracy better than 8% in the x and 
Q range studied for both CTEQ4M and CTEQ4LQ (see Fig. This is very good as far 
as DVCS studies are concerned since Eq. (|T0|) is an all order statement and thus one can 
use the NLO evolution of diagonal gluon densities to make NLO predictions for DVCS! As 
mentioned before, we have chosen X2 = Xi — A = 10~® for xi = 10~^ and X2 = 10"'' for 
Xi = 10~^. As one can see there is no big difference in the ratio of the densities for the two 
A values studied which is in agreement with Radyushkin's statement that Eq. (|l^) should 
hold as long as A is much smaller than 1. 

Finally, we would like, in light of the previous findings, comment on the NLO evolution 
in the DGLAP region. Given the fact that Eq. (p!0| , pJ]) was based on general arguments in 



23| and that it holds in LO evolution with a certain degree of accuracy, we predict that the 
NLO evolution will not change the above relations, in other words, that the NLO evolution 
of the nondiagonal gluon distribution can be predicted to a similar accuracy by the NLO 
evolution of the diagonal gluon distribution. We base this statement on the results of the 



^The same behaviour was observed by Golec-Biernat [25 
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above analysis and the observations of that the NLO corrections of the nondiagonal 
evolution should be in the same direction as in the diagonal case, which reduces the LO 
results, and of the same magnitude. The former statement is due to the observations made 
in 1^ that if, in the nondiagonal case, the NLO corrections were in the opposite direction, 
which would lead to a marked deviation from the LO results, compared to the diagonal case, 
the overall sign of the NLO nondiagonal kernels would have to change for some A 7^ since 
in the limit A ^ we have to recover the diagonal case. This occurance is not likely for 
the following reason: First, the Feynaman diagrams involved in the calculation of the NLO 
nondiagonal kernels are the same as in the diagonal case, except for the different kinematics, 
therefore, we have a very good idea about the type of terms appearing in the kernels, namely 
polynomials, logs and terms in need of regularization such as \n(z) x . Secondly, the 

kernels, as stated before, have to reduce to the diagonal case in the limit of vanishing A 
which fixes the sign of most terms in the kernel, thus the only type of terms which are 
allowed and could change the overall sign of the kernel are of the form 

-/(xi/yi,A/yO (13) 

yi 

which will be numerically small unless ?/i ~ A in the convolution integral of the evolution 
equations. Moreover, we know that in this limit the contribution of the regularized terms 
in the kernel gives the largest contributions in the convolution integral and therefore sign 
changing contributions in the nondiagonal case would have to originate from regularized 
terms. This in turn disallows a term like Eq. |13| due to the fact that regularized terms are 
not allowed to vanish in the diagonal limit, since the regularized terms arise from the same 
Feynman diagrams in the both diagonal and nondiagonal case. Therefore, the overall sign 
of the contribution of the NLO nondiagonal kernels will be of the same as in the diagonal 
case . In addition, the magnitude of the correction should be of the same magnitude as in 
the diagonal case since one has the conditions R{Q) > 1 P J20|J23[] , R{Q) < 1.5 for the 
gluon distribution and the fact that the LO results at high Q are already fairly close to the 
upper bound. This forces the NLO corrections in the nondiagonal case not to exceed the 
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diagonal corrections by a factor of 3 or so, lest it violates the boundary conditions for R{Q). 



V. CONCLUSIONS 



In the above, we examined predictions made in Ref. about relationships between 
nondiagonal and diagonal parton distributions in the DGLAP region based on certain models 
for nondiagonal parton distributions in the normalization point. We found that the evolution 



does not destroy the validity of Eqs. (|10|,|TT]) for all x and Q > Qo for both CTEQ4M and 
CTEQ4LQ in the case of gluons and for xi > 2^3A for quarks. Therefore, we conclude that 
Eqs. ([TO|JTT|) do supply a reliable approximation of nondiagonal parton densities for A ^ 1. 
We also conclude that we can hardly see a variance of the accuracy of predictions (|lO| , [IlD 
for different initial distributions. The accuracy is slightly better for the one which supplies 
a less steep gluon density at small x. 

Based on these results and the results from Ref. [0], we predicted the NLO evolution of 
the nondiagonal gluon distribution to be within ^ 20% of the diagonal gluon distribution 
for the above made ansatz and for a large range of Q. 
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After this work was complete, we learned that Miiller et al. |2^ had performed a numer 



ical study of the NLO effects for the non-singlet and singlet distributions and found that 
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the corrections were within 10 — 30% percent of the LO result confirming our statements on 
the NLO corrections. 
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FIG. 2. R{Q) is ploted as a function of Q at A = 10 ^ and A = 10 ^. The input distributions 
are CTEQ4M and CTEQ4LQ. 
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FIG. 3. c/(xi, A,Q)/(xi - A/2)G(a;i - A/2,Q) as a function of xi for Q=1.6, 2.5, 5.3, 10.0 

GeV for A = IQ-"^ and A = 10~^. The initial distribution is CTEQ4M. 
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FIG. 4. c/(xi, A,Q)/(xi - A/2)G(a;i - A/2,Q) as a function of xi for Q=0.7, 2.8, 5.9, 10.0 

GeV for A = IQ-"^ and A = IQ-^. The initial distribution is CTEQ4LQ. 
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FIG. 5. q{xi,A,Q)/{xi - A/2)Q(xi - A/2,Q) as a function of xi for Q=1.6, 2.5, 5.3, 10.0 

GeV for A = IQ-"^ and A = IQ-^. The initial distribution is CTEQ4M. 
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FIG. 6. q{xi,A,Q)/{xi - A/2)Q(xi - A/2,Q) as a function of xi for Q=0.7, 2.8, 5.9, 10.0 

GeV for A = 10"^ and A = IQ-^. The initial distribution is CTEQ4LQ. 
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